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Abstract 

This is a concise review of the complex, real and quaternion real Ginibre random 
matrix ensembles and their elliptic deformations. Eigenvalue correlations are 
exactly reduced to two-point kernels and discussed in the strongly and weakly 
non-Hermitian limits of large matrix size. 

• i— i 
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^ . 18.1 Introduction 

The study of eigenvalue statistics in the complex plane was initiated in 1965 
by Ginibre |Gin65| more than 40 years ago who introduced a three fold family 
of Gaussian random matrices (complex, real and quaternion real) as a math- 
ematical extension of Hermitian random matrix theory. Although no physical 
applications of the theory were in sight at that time, Ginibre expressed the 
hope that 'the methods and results will provide further insight in the cases of 
physical interest or suggest as yet lacking applications'. Nowadays, statistics of 
complex eigenvalues have many interesting applications in modeling of a wide 
range of physical phenomena. They appeared in the studies of quantum chro- 
modynamics (see Chap. 32), dissipative quantum maps [Gro88| and scattering 
in chaotic quantum systems (see Chap. 34), growth processes (see Chap. 38), 
fractional quantum-Hall effect [DiF94] and Coulomb plasma [For97], stability 
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of complex biological |May72| and neural networks |Som88j . directed quantum 
chaos in randomly pinned superconducting vortices |Efe97| . delayed time se- 
ries in financial markets [Kwa06], random operations in quantum information 
theory [Bru09], and others. 

This Chapter gives an overview of the three Ginibre ensembles and their 
elliptic deformations. We tried to keep our exposition self-contained providing 
hints of derivations. Some of the derivations included are new. We did not have 
space to cover the chiral extensions of the Ginibre ensembles, only mentioning 
briefly the chiral companion of the complex Ginibre ensemble [Osb04| . This 
topic is partly covered in Chap. 32. The real and quaternion real (qu-r) com- 
panions were solved only recently, see |Ake09c| and [Ake 05aj . There are also 
important non-Hermitian ensembles of random matrices relevant in the context 
of quantum chaotic scattering which are only mentioned here (but see Chap. 
34 for a summary of results and the survey paper |Fyo03| for details). Also, we 
will not discuss the complete classification of non-Hermitian matrix ensembles 
depending on the action of a few number of involutions [Ber021 Mag08| . 

On a macroscopic scale, all three Ginibre ensembles exhibit similar patterns 
of behavior with a uniform distribution of eigenvalues and sharp (Gaussian) fall 
in the eigenvalue density when one transverses the boundary of the eigenvalue 
support. The similarities extend to the microscopic scale as well but only 
away from the real line where all three ensembles exhibit a cubic repulsion of 
eigenvalues. In the vicinity of real line and on the real line their behavior is 
very different due to the differences in symmetries. The eigenvalue correlation 
functions have either determinantal (complex) or pfaffian (real and qu-r) form 
with the kernel being almost identical far away from the real line, again the 
differences coming from a pre-exponential factor describing the transition from 
the real line into the bulk of the spectrum. 

One expects the eigenvalue statistics provided by the Ginibre ensembles to 
be universal within their symmetry classes. Establishing such universality is 
an open and challenging problem, especially for the real and qu-r ensembles. 
There is some evidence for universality in the complex case where every solved 
model led to Ginibre's form of correlations, including complex normal matrices 
where the universality of Ginibre's correlations was proved in a general class of 
matrix distributions [Ame08j. This result can be applied to complex matrices 
(118. 2. 4p as for this class of ensembles the induced eigenvalue distribution does 
not differ from the one for normal matrices [Oas97j . see also Chap. 38. 



18.2 Complex Ginibre Ensemble 



Measure, change of variables The complex Ginibre ensemble is defined 
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on the space of complex N x N matrices by the probability measure 

d/i(J) = exp(-TrJJ t )|DJ|. (18.2.1) 

Here DJ = nfj=i(d<A?dJjjV27r) is the (exterior) product of the one-forms in 
matrix entries and |DJ| is the corresponding Cartesian volume element. With 
probability one, the matrix J has N distinct eigenvalues Zj. On ordering the 
eigenvalues in an arbitrary but fixed way, one can think of them as of random 
variables. The corresponding joint probability distribution function (jpdf) can 
be obtained by changing variables in (|18.2.ip . This can be done in several ways, 
we outline here a calculation due to Dyson |Meh04] . 

On making use of the Schur decomposition, J can be brought to triangular 
form J = f7(A + A)U~ 1 . Here U is unitary, A = diag(^i, . . . ,zn) and A is 
strictly upper-triangular. It is apparent that U can be restricted to the space 
of right cosets \J (N) /XJ (1) N . The variations in J are related to those in A, A 
and U by the formula U~ l dJU = dA + dA + dM, where 

(dM)y = (U^dU^Zj - Zi ) + J2( U ~ ldU )ik&kj -J2 A il( U ~ ldU h ■ 

k<j l>i 

The volume form DJ does not change on conjugation by unitary matrices. 
Hence DJ = D(A + A + M) and one gets the Jacobian of the coordinate 
transformation by multiplying through the matrix entries of dA + dA + dM 
using the calculus of alternating differential forms. This leads to the important 
relation 

d /U (J) = CVe- Tr(AAt+AAt) Y[ |z; - z/ |D[/| |DA| |DA| , (18.2.2) 

l<i<j <N 

where DC/ = f\i<j (U~ 1 dU)ij(U~ 1 dU)ji is a volume form on the coset space. 
The density function on the rhs is symmetric in the eigenvalues of J. It does not 
depend on U and is Gaussian in A, and these two variables can be easily inte- 
grated out. Thus, if f(zi, . . . , zn) is symmetric in eigenvalues of J then [Gin65j 
/ fdfjt = J d 2 z x • • •/ d 2 z N P N (z 1 , . . .,z N )f(zi, ...,z N ) where d 2 z = \dzdz*\/2 
is the element of area in the complex plane and 

P(z 1 ,...,z N )= * - e -E£i'**' a J] \zi-Zjl 2 (18.2.3) 

n llj=o3- l<i<j <N 

is the symmetrized jpdf of the eigenvalues of J. 

In contrast to the Hermitian matrices, the above calculation of the jpdf is 
not easily extended to other invariant distributions. E.g., it breaks down if Jjt 
in (118.2.ip is replaced by a higher order polynomial in Jjt [Fei97] as U, A and 
A do not decouple in that case, although it still works well for (see Chap. 38) 



d/i( J) tx exp ( - Tr JJ f - Re Tr <3?(J)) |D J 



(18.2.4) 
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with ^(J) being a potential that ensures existence of the normalisation integral. 
We will discuss in detail the special case of (|18.2.4|) with $(J) = J 2 later. 

Correlation functions and orthogonal polynomials The eigenvalue 
correlation functions are marginals of the symmetrized jpdf, 

Rn(zi, ...,z n )= N ' [ d 2 z n+1 ... [ d 2 z N P(z 1 , ...,z N ), (18.2.5) 



(N -n 

normalized to / d 2 z 1 ... j d 2 z n R n (zi, ...,z n ) = N(N - 1)...(N - n + 1). The 
one-point correlation function is just the density of eigenvalues p(z) = £^ • 5^ (z— 
Zj) averaged over the ensemble distribution, Ri(z) = (p(z)), so that if no is 
the number of eigenvalues in D then (no) = J D d 2 zRi(z). 

The eigenvalue correlation functions for the complex Ginibre ensemble can 
be found in a closed form. The corresponding calculation is almost identical 
to that for Hermitian ensembles. For the purpose of future reference we shall 
outline it in a slightly more general setting. 

Let Pm(z) be the monic orthogonal polynomials associated with weight func- 
tion w(z) > in the complex plane, i.e., J d 2 z w(z) p m (z)p n (z)* = h n 5 m ^ n and 
Pm(z) = z m + . . .. Then 

1 N 

P(*i, ■■-,zn)= A7i tt N—i 7 II II N - %| 2 (18-2-6) 
^ ! llz=0 h l j=l i<j 

is a probability density in symmetric with respect to permutations of Zj. Re- 
calling the Vandermonde determinant T\i<i<j<N( z i ~ z j) = ( ^ e ^( z i f ~'')i<i,j<N = 
det(pN-j(zi))±<ij<N, one obtains the important determinantal representation 
of the jpdf: 

P(zi, ...,z N ) = ^ det(K(z u Zj ))% =1 (18.2.7) 

with 

N-l ( \ ( \* 

K( Zl ,z 2 ) = v^v^) £ MZl f" (Z2) . (18.2.8) 

n=0 

The kernel K is Hermitian, K(z\,z 2 ) = K(z2, Z\)* , and J d 2 z K(z,z) = N, 
J d 2 z K {z\, z)K(z, z 2 ) = K{z\,z 2 ). Hence by Mehta's 'integrating out' lemma 
[Meh04] (see also the relevant section in Chap. 4) 

R n (z 1 ,...,z n ) = det(K{z i ,z j ))l J=1 . (18.2.9) 

In particular, the one- and two-point correlation functions are given by R\{z) = 
K(z, z) and R 2 (zt, z 2 ) = K(z 1 , z 1 )K(z 2 , z 2 ) - \K(z 1 ,z 2 )\ 2 . 

In the complex Ginibre ensemble w(z) = e - ' 2 ' . In view of the rotational 
symmetry the power functions are orthogonal. Thus p n (z) = z n , h n = im\ and 

K{zijZ2 ) = I e-H^I 2 -il^l 2 y ifi^L. (18.2.10) 

n=0 
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The sum on the rhs is the truncated exponential series and can be expressed in 
terms of the incomplete Gamma function: X^o* 7T = &X ^0^' X )/T(N). The 
saddle-point integral T(N,x) = dt e - *^ -1 comes in handy for asymptotic 
analysis. We have 

R 1 (z) = - T{ ^ l f } ~ le(y/N-\z\), iV-oo, (18.2.11) 

7T 1 [N ) 7T 

where is the Heaviside function, meaning that on average most of the eigenval- 
ues are distributed within the disk \z\ < V~N with constant density in agreement 
with the Circular Law |Gir851 ITao08j . The expected number of eigenvalues out- 
side this disk ~ \f^- The eigenvalue density in the transitional region around 
the circular boundary, as found from the integral T(N,x), is given in terms of 
the complementary error function 

Rt{(VN + x) e*) ~ — erfc(V2z), as N -> cxd, (18.2.12) 

|For991lKan05a] . Recalling that erfc(x) ~ 2 0(— x) + e~ x2 / (t/ttx) for large real 
\x\, one concludes that the density vanishes at a Gaussian rate at the edge. 

The kernel in (|18.2.10p has a well defined limit as N — » oo and |zi 2I = 0(1), 
leading to a simple expression for the correlation functions in this limit [Gin65j 

R n (zi, . . . , z n ) ~ -L e" ^ ^ det(e**?)? J=1 . (18.2.13) 

Equation (|18.2.13p describes eigenvalue correlations at the origin zq = and on 
the (local) scale when the mean separation between eigenvalues is of the order 
of unity. It also holds true (locally) at any other reference point inside the 
disk \z\ < y/~N, see, e.g. |Bor09| . In particular, Rzizii&z) — -^(l — e - ' 21-22 ' ). 
Note that the dependence of R n on the reference point disappears in the limit 
N — > 00 and complete homogeneity arises. The eigenvalue correlation functions 
at the edge of the eigenvalue support can also be found, see |For991 lBor09j . 

Gap probability and nearest neighbor distance Consider the disk 
D of radius s centered at zq and denote by xd its characteristic function. Define 
H(s;zq) to be the conditional probability that given one eigenvalue lies at zq 
all others are outside D. With xd being the characteristic function of D, 

H (s;z ) = — — — / d 2 z 2 --- / d 2 z N P(z ,z 2 , - ■ ■ ,z N )Y\(l - XD(zk)) ■ 
Ri{z ) J J * = A 2 

(18.2.14) 

When s gets infinitesimal increment 6s, the decrement in H, H(s, Zq) — H(s + 
6s, Zo) , is the probability for the distance between the eigenvalue at zq and 
its nearest neighbor to lie in the interval (s, s + 5s). Therefore p{s;z$) = 
— ■^H(s;zq) is the density of nearest neighbor distances at Zq. 
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For the complex Ginibre ensemble H(s; 0) can be easily computed with 
the help of the Vandermonde determinant and Andreief-de Bruijn integration 
formula (Eq. (8.3.41) in Chap. 8). The rotational invariance of the weight 
function e - ' 2 ' ensures that the monomials z m stay orthogonal when integrated 
over the disk D, leading to |Gro88] 

N-i / 



T(n+ l,s 

, r(« + i) 

n=l v ' n=l 




The product on the rhs converges quite rapidly in the limit N — > oo and H (s; zq) 
has a well defined limit. In this limit H(s, 0) = l-s 4 /2+s 6 /6-s 8 /24+0(s 10 ) for 
small s and [Gro88| H(s,0) = exp[— s 4 /4 — s 2 (lns + O(l))] for large s. Taking 
the derivative, one gets the cubic law of the eigenvalue repulsion: p(s;0) = 
2s 3 + 0(s 5 ) for small s. 

The small-s behavior of H(s, zq) can also be obtained by expanding the 
product in (118.2.140 . E.g., on retaining the first two terms H(s,zq) = 1 — 
Ri(z ) f\z-z \<s ^ Z Rz( z 0' z ) + 0(s A ). This approach is general and does not 
rely on the rotational invariance. The universality of Ginibre's correlations 
then implies universality of the cubic law of the eigenvalue repulsion, which 
can also be verified directly for complex ensembles with known jpdf [Oas97| . 
Retaining all terms leads to a Fredholm determinant, see Chap. 4, giving access 
to various gap probabilities beyond the Ginibre ensemble [A~ke09bJ . 



18.3 Random contractions 



Let U be a unitary matrix taken at random from the unitary group U(M + L). 
Denote by J its top left corner of size M x M. The matrix J is a random 
contraction: with probability one, it has all of its eigenvalues inside the unit disk 
\z\ < 1. The jpdf of eigenvalues of J was computed in |ZycOO| and is given by 
(|18.2.6p with wl(z) = (1— |2:| 2 ) L_1 G(1 — |^|). Since wl{z) is rotation invariant, 
the associated orthogonal polynomials are again powers p m (z) = z m and the 
normalization constants are easily computed in terms of the Beta function, 
h m = TrB(m + 1,L). Applying the orthogonal polynomials formalism, one 
obtains the correlation functions in the determinantal form (|18.2.9p with 

T M-l 

K{z 1 ,z 2 ) = -{l-\z 1 \^{l-\z 2 \^ V 

7T 

The truncated binomial series for (1 — x)~( L+1 ^ on the rhs can be expressed in 
terms of the incomplete Beta function I x {a, b) = B ^ a b ^ Jq i a-1 (l — t) b ~ 1 dt via 

the relation I X (M,L + 1) = 1 - (1 - x) L+1 *Em=o ^m") ^ This leads to a 



(L + m\ 
\ m J 
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useful representation 



L (l-biP)^(l-hl 2 )^ 



K(zu z 2 ) = - i li - ZlZ *)Lii i 1 " ( M ' L + !) J • ( 18 - 3 - 16 ) 

There are several asymptotic regimes to be considered in the context of M x M 
truncations of random unitary matrices of an increasing dimension M + Z0. 
The simplest one is the limit when M stays finite and L — > oo. In this regime 
the jpdf becomes Gaussian and one immediately recovers the complex Ginibre 
ensemble. In fact one can allow M to grow with L and still recover the Ginibre 
ensemble provided that M <C L; see .)ia()(i| and references therein for bounds 
on the rate of growth of M for the Gaussian Law to hold. In the opposite 
case when the size of the truncation increases at the same rate as the overall 
dimension M + L, one has two distinct regimes: (i) M, L — * oo, M/L = a > 0, 
and (ii) M — > oo, L is finite. 

We start with (i) which is the limit of strong non-unitarity. In this limit the 
eigenvalues are distributed inside the disk \z\ 2 < with density |Zyc00| 

, > _ L 1-I N 2(M,L + 1) M 1 ^ ( a |2 
llZj "vr (l-k| 2 ) 2 ~ Tra (1 - |z| 2 ) 2 Vl + a |Z| 



in the bulk and Ri[Jj^ + f=) - exfc (V5 ^ a?) at the edge. 

Modulo a simple rescaling, the edge density is the same as in the complex 
Ginibre ensemble. It can also be seen that the average number of eigenvalues 
outside the boundary ~ y / M(l + a)/(2ii). 

One can also find the eigenvalue correlation functions in the bulk. At the 
origin this task is especially simple. Scaling z by V ttL, which is the mean 
distance between the eigenvalues at the origin, and noticing that for u in a 
bounded region in the complex plane ^m=o flT) (x)™ — e " ™ s t ron g 
non-unitarity limit, one concludes that the rescaled correlation functions are 
given by Ginibre's expression f|18.2. 13 j) . Although the eigenvalue distribution 
now is not homogeneous, after appropriate rescaling (|18,2.13j) describes the 
eigenvalue correlations at any point in the bulk |Zyc00| . 

The gap probability at the origin, H(s;0), can be computed by exploit- 
ing the rotational symmetry, in the same way as for the Ginibre ensemble, 
H(s;0) = n m =i (l — I s 2(m + l,L)). After appropriate rescaling one obtains 
the same expression as in the Ginibre ensemble (cf. (| 18.2. 15 j) ) . H ^-^;0^ ~ 



lr The integral I x (a,b) is convenient for finding the relevant limits by the Laplace method. 
E.g., one finds that I x (a,b) ~ Q(x — ^rj) as a,ii -> oo. The remainder term here is ex- 
ponentially small when x is away the transitional point xq = ^ " a \ ■ In the transitional 

region I XQ+ _^(a,b) ~ 1 — | erfc ( ^ 1+ J^ 1 t) ■ Another asymptotic relation of interest is 



I 1 _v(a,b) ~ r(6, y)/T(y) which holds for positive y in the limit a — + oo, b is fixed. 
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Um=i T ( m + !) s 2 )/T(m + 1) in the limit M, L oo, M/L = a, leading to the 
cubic law of of eigenvalue repulsion in the bulk. 

Now, consider the limit of weak non-unitarity M — > oo, L is finite. In this 
limit the eigenvalues of J lie close to the unit circle, with the magnitude of the 
typical deviation being of the order 1/M. Scaling z accordingly, one finds the 
eigenvalue density in the limit of weak non-unitarity |ZycOO| : 

Ri(l~T7^) - — fy"> L1 f e-W-dt, M — oo and L is fixed. 
V M J 7T (L-1)\J Q 

Setting Zj — (l — f|) e 1 ^ 0-1 " 1 "^, one finds the correlations in this limit: 



Interestingly, a different random matrix ensemble, J = H + vyW, leads to the 
same form of the correlation functions (after appropriate rescaling) as on the 
rhs above |Fyo99|. Here H is drawn from the GUE, 7 > and W is a diagonal 
matrix with L l's and M zeros on the diagonal, with M ^> 1 and finite L. The 
above equation is a particular case of a universal formula describing correlations 
in more general ensembles of random contractions and non-Hermitian finite rank 
deviations from the GUE, see |Fyo03| and Chap. 18 for discussion and results. 



18.4 Complex elliptic ensemble 



Let H\, H2 be two independent samples from distribution dfi(H) = e \DH\ 
on the space of Hermitian N x N matrices. Then J = y/T~+ tH\ + — tHi 
is a random matrix ensemble interpolating between the (circular) Ginibre en- 
semble (r = 0) and the GUE (r = 1). We only consider the interval < r < 1. 
The matrix J is complex Gaussian, 

dfi(J) oc exp { Tr TjJ f + -(J 2 + J t2 )l \ |DJ| . (18.4.17) 

I 1 — t z L 2 J J 

The jpdf, P(zi, . . . , zn), can be obtained by bringing J to triangular form, 
as in Section 118.21 The resulting expression is given by (|18.2.6p with weight 
function w T (z) = 2 exp{ — + ^jjfsy^ } ■ The associated orthogonal 

polynomials are scaled Hermite polynomials [DiF94]. Indeed, by making use 
of the integral representation H n (z) = ^e 2zi_ ' f~( n+1 ) dt, for the Hermite 
polynomials H n (z) one can easily verify that 



18.4. COMPLEX ELLIPTIC ENSEMBLE 



9 



The monic polynomials p n (z) are easily found, p n {z) = (r /2) n / 2 H n {z / 'y/2z) , 
h n = n\, and on applying the orthogonal polynomial formalism, one obtains 
the correlation functions in the determinantal form (|18.2.9p with 

K(z 1 ,z 2 ) = w 1 T / 2 (z 1 )w 1 J 2 (4) J2^H n Uy)H n (^=) . (18.4.19) 

n=0 ' 

The sum on the right is the truncated exponential series in Mehler's formula 

A quick comparison of (j!8.4.19j) with Mehler's formula convinces that the den- 
sity K(z, z) is constant in the limit A — > oo. More care is needed to determine 
the boundary of the eigenvalue support. To this end, another integral for the 
Hermite polynomials comes in handy, H n (z) = (±2i) n J^°° e~(* ±lz ) t n dt/y/n. 
On substituting this into (|18.4.19p and after a simple change of variables, one 
writes the kernel in a form suitable for asymptotic analysis 

1 f |zil 2 + IZ2I 2 — 2z\Zn ~| 

K(zu^) = 2(1 -r') J X (18A20) 

> £> r ( " )/+W rWAr r T ( ( ;,"" 2)> ' (1M ' 21) 



with/±( 9 ) = v ^pexp{-(,ViV(T±7)-^gi ! ) }■ 

For large with 1 — r > uniformly in N and x,y = 0(yN), the functions 
can formally be replaced by delta-functions, f + (v) = S(y — — ) and 

f~(u) = S(u + ^ — L=J. Such a replacement can be justified by deforming 

the u-integral into the complex plane to pick up the sharp peak of f + (v) along 
the imaginary axis. This gives a simpler expression for the eigenvalue density 
Ri{z) ~ (tt(1 - t 2 ))~ 1 T(N, jt ^ 7 + T ^- ? )/T{N), cf (|18.2.11|) . Hence, for 
large A the eigenvalues density is l/vr(l — r 2 ) inside the ellipse with half-axes 
v / A(l ± t) along the x and ?/ directions, in agreement with Girko's elliptic law 
[Gir86]. It falls to zero exponentially fast when one transverses the boundary 
of the ellipse. 

For z\, Z2 inside the ellipse and such that \z\ — z<i\ = 0(1) the integral in 
(I18.4.2ip converges to 1 as A — > oo, and after a trivial rescaling one obtains 
the same expression (I18.2.13P for the correlations as in the circular case (at the 
origin this readily follows from Mehler's formula). 

The limit A — > oo, 1 — r > is the limit of strong non-Hermiticity. There 
is another important limit, A — > oo and 1 — r = a 2 /A, the so-called limit of 
weak non-Hermiticity |Fyo97| describing the cross-over from Wigner-Dyson to 
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Ginibre eigenvalue statistics. In this limit the eigenvalues of J can be thought of 
as those of \/2H\ (the Hermitian part of J) displaced from the real axis by the 
"perturbation" term \-^=H2 (the skew-Hermitian part of J). The eigenvalues 

xj of V2Ht fill the interval (-2-/ZV, 2yfN) with density p sc (x) = ±s/N - x 2 /A 
(Wigner's Semicircle Law). A simple perturbation theory calculation Fyo97 
gives the density of eigenvalues z = x + \y of J in the factorized form p sc (x)p(y), 
with the Gaussian distribution p(y) = y N/2ira 2 exp(— Ny 2 /2a 2 ) of the dis- 
placements. For such a calculation to be well defined, the width of this dis- 
tribution should be much smaller than the (mean) eigenvalue spacing l/p sc (x) 
of the unperturbed eigenvalues, making it natural introducing the control pa- 
rameter a(x) = ap sc (x)/y/~N. Formulas (j!8. 4.20 j )- (|18. 4.21 j) make it possible 
to go beyond the perturbation theory and compute the eigenvalue density and 
correlations in the limit of weak non-Hermiticity exactly. In this limit the func- 
tion f + (v) is singular and the same as in the limit of strong non-Hermiticity, 
however the function f~(u) is not. Let Zj = x + Q/p sc (x). Then, to the lead- 
ing order, f~(u) ~ ix~ l l 2 a exp (ua — ^ 7§j ) 2 anc ^ ^ ae integral in (|18.4.2ip 

~ 2^pKa exp{ ^ 1 4a i 2 ' > } J*q e~ 7T ' 2a2u2 cos(iru((i — CI)) du, where a = a(x) is the 
control parameter defined above. One then obtains the scaled correlation func- 
tions Rn(Ci, . . . , Cn) = (p 2 sc (x))~ n R n (zi, ...,z n ) in the form |Fyo97| : 

Mb, ■ ■ ■ , Cn) * (~j=A ^ (Im ^ )2 det (J^'^ 2 cos(7r«(C< - CJ)) d«) . 

(18.4.22) 

It is easy to check that R n interpolates between the Wigner-Dyson correlations 
(a — > 0) and Ginibre's (o — > oo), and the above equation allows one to study 
the cross-over from one set of eigenvalue statistics to the other, see |Fyo98| . 
The eigenvalue density is 

1 (ImQ 2 f 1 2 2 2 

i?i(C)^^^e ^~ I e~ n a u cosh(27ruImC) du , (18.4.23) 

correcting the perturbation theory result. Interestingly, the distribution of the 
scaled imaginary parts of eigenvalues as given in (I18.4.23j) appears to be uni- 
versal. A supersymmetry calculation shows |Fyo98| that it does not depend 
on the details of the Hermitian and skew-Hermitian parts of J and extends to 
non- invariant matrix distributions (like the semicircular and circular laws). It 
is an open challenging problem to find a proof of this universality satisfying the 
rigor of pure mathematics. Staying within the class of invariant distributions, 
the eigenvalue density and correlations appear to be also universal: it was ar- 
gued in [Ake03j that the weakly non-Hermitian limit in the ensemble (|18.2.4[) 
at the origin is also described by (j!8.4.22|) . 

In conclusion, we would like to mention briefly two topics related to the 
complex elliptic ensemble. One is the recent studies of edge scaling limits 
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|Gar02l IBen09j . By scaling r with N so that 1 - r = a/A 1 / 3 |Ben09 ] one 
gets access to the crossover from Airy (aCl, GUE) to Poisson edge statistics 
(a 3> 1, Ginibre). And the other is the chiral extension J = ( iyl t+ M £t 
of the complex elliptic ensemble. Here A and B are two independent samples 
from the Gaussian measure with density e - TrA '' A Q n the space of complex 
(JV + v) x N matrices, i/>0, 0</i<l. This ensemble can be studied along 
the same lines as (I18.4.17p . The computation of the jpdf of eigenvalues [Osb04] 
is though more involved and leads to (|18.2.6p with weight function w(z) = 

\z\ 2u+2 e 4 '' 2 ^ z +z ^ Ku C^ \ z \ 2 ) where K u is a modified Bessel function with 
the associate orthogonal polynomials being scaled Laguerre polynomials [Osb04l 
IAke05a| . At // = we recover the Wishart (Laguerre) ensemble of Hermitian 
matrices and the corresponding weakly non-Hermitian limit N — > oo and a = 
N/j, 2 = O(l) describes its neighborhood. The eigenvalue correlations in this 
limit were computed in [Osb04], with the answer being somewhat different 
from (|18.4.22p . This is not surprising given that the two ensembles belong to 
different symmetry classes. By letting a — > oo one obtains the correlations in 
the chiral ensemble in the regime of strong non-Hermiticity, see, e.g., [Ak e05b| . 

18.5 Real and quaternion-real Ginibre ensembles 

Measure, change of variables Restricting to even dimensions N, we will 
treat the real and quaternion real (qu-r) ensembles in a unifying way. Consider 
the normalized measure 

d/i(J) = exp(-- Tr J.P) |DJ| (18.5.24) 



on the space of N x N matrices subject to the constraints 

J T = in the real case 

ZJ T Z T = in the quaternion-real (qu-r) case . 



(18.5.25) 



where Z = (-1 o) ^ s ^ ne symplectic unit. Here DJ = nfj=i( c ^ij/v / 2vr) 

is the (exterior) product of the one- forms in matrix entries and |DJ| is the 
corresponding Cartesian volume element. 

As follows from (|18.5.25p . both real and qu-r matrices have non-real eigen- 
values occurring in pairs z and z* , and any real eigenvalue of a qu-r matrix 
has multiplicity > 2. Hence, the probability for a qu-r matrix drawn from the 
distribution (|18.5.24|) to have a real eigenvalue is zero. 

As with complex matrices, the jpdf of eigenvalues can be obtained from a 
Schur decomposition. Because of the symmetries, it is convenient to work with 
matrices partitioned into 2x2 blocks. Ignoring multiple eigenvalues, we can 
bring J to block-triangular form, J = U(A + A)C/ _1 where A is block-diagonal 
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and A has nonzero blocks only above A. The matrix U is orthogonal in the 
real case and unitary symplectic in the qu-r case and can be restricted to the 
space of right cosets, 0(N)/0(2) N / 2 and USp(iV)/USp(2) 7V / 2 respectively. 
The variations in J are related to those in U, A, A by 

U~ l dJ U = dA + dA + [U~ x dU, A] + [U~ l dU, A] , (18.5.26) 

with [A, B] = AB — BA. The matrix U~ 1 dU is skew-symmetric in the real case, 
{U^dU f = -(U- x dU), and anti self-dual in the qu-r case, Z(U' 1 dU) T Z T = 
-OJ- l dU). 

To find the Jacobian associated with the change of variables from J to 
U, A, A, we multiply entries of the matrix on the rhs in (|18.5.26p using the 
calculus of alternating differential forms. The matrices dA and dA yield the 
volume forms for A and A, respectively. Due to the triangular structure, the 
matrix [U~ 1 dU, A] does not contribute, and the matrix [U~ 1 dU, A] yields the 
coset volume form times a factor depending on the eigenvalues Xj of A, 

^(U^dUA - MJ^dU)^ = [jV _1 d^)ii(Aj - Ai) , 

with the dashed product running over non-zero entries in the lower triangle of 
U~ 1 dU. On gathering all terms, we arrive at 

d/i(J) = e -| Tr ( AAt+AAt ) |DA| |DA| inV -1 dEOy(Aj - (18.5.27) 

with DA = W (dAjj/\/27r) and DA = n (dA^/ \/2n), the products running over 
nonzero entries. The jpdf of eigenvalues follows from (|18.5.27p on integrating 
out all auxiliary variables. Since the density function does not depend on U, the 
{/-integral gives the volume of the coset space. The integral in A is Gaussian, 
and this variable can be integrated out with ease as well. Thus we are left with 
the problem of integrating over the N/2 (2 x 2)-blocks appearing in A keeping 
the eigenvalues in each block fixed. These are precisely the eigenvalues of the 
matrix J and our problem reduces to 2 x 2 matrices. 

Dimension N = 2 and the jpdf The generic form of a qu-r 2x2 
matrix is J = ( __^* a b * ) with a and b complex. The eigenvalues are complex 
conjugate, Ai^ = Re a ± iy 7 (Im a) 2 + \b\ 2 , and we choose ImAi = — ImA2 > 0. 
Thus, Jdfi(J) 5 2 (A X - Re a - iv^Ima) 2 + |6| 2 ) = ^M/9(A X ) 2 with /<?(A) as 
in (|18.5.29p . The rhs is normalized wrt integrating over Ai in the half-plane 
ImAi > 0. 

Turning to real 2x2 matrices J = (™ ^) with a, b, c, d real, the eigenvalues 
are Ai^ = \{a + d ± \/ (a — b) 2 — 46c), and we now have two possibilities (i) 
Ai = A2* (choosing ImAi = — ImA2 > 0), or (ii) both Ai and A2 are real 
(choosing Ai > A2). In the first case we obtain [Leh911 ISoinT)8] / dn(J)5 2 (X 1 - 
\[a + d - i^/Abc- {a-b) 2 }) = ^M/ r (Ai) 2 with / r (A) as in (118.5.291) . In 
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the 2nd case (Ax,2 real) the corresponding average over the product of two 
5-constraints results in %^2l f (A x ) f r (X 2 ) |Som08j . 



2V2tt 

Consequently, the joint distribution of eigenvalues for both ensembles real 
and qu-r can be written in the form |Gin651 IMeh044 ISom08] 

N 

d f x(\ 1 ,X 2 ,...,X N ) = C N H (X{ - Xj) JJ /(Aj) dAidAa • • • dAjv , (18.5.28) 

l<i<j<N i=l 

where 

(f r (X) 2 = erfc ( e"5( A2+A * 2 ) in the real case 
/(A) 2 = I K V2 J (18.5.29) 

(fi(X) 2 = | A - A* | e"l A l 2 in the qu-r case, 

with /(A) = /(A*) > in both cases. We have to put the eigenvalues in such an 
order that d/x > 0. Thus we consider each case with m complex conjugate pairs 
of eigenvalues (m = 0, 1, 2, N/2) separately, arranging the Aj's in (|18.5.28p 
in the following order: for m = (all real) Ai > A2 > ••• > An, for m = 1 (one 
complex conjugate pair Ai = A2*) ImAi > Im A2, A3 > A4 > ... > Ajv, for m = 2 
(two complex conjugate pairs Ai = A2* and A3 = A|) ImAi > ImA2,ReA2 > 
Re A3, ImA3 > ImA4, A5 > ... > Xn, etc.. Summing and integrating over 
all cases and ranges yields total probability 1. In the qu-r case we only have 
m = N/2 complex conjugate pairs, and in the real case the probabilities of 
finding m complex conjugate pairs have been calculated in |Kan05b| IFor07] . 

The jpdf (|18.5.28p shows immediately the repulsion behavior of the eigen- 
values. Due to the Vandermonde it contains for two eigenvalues in the upper 
halfplane the factor |Ai — X 2 \ 2 which means cubic repulsion in the distance, the 
additional power coming from the two-dimensional volume element [Gro88j. 
For two eigenvalues exactly on the real axis one has the factor |Ai — A2I and 
thus linear repulsion like in GOE. This applies only in the real case. In this case 
due to the factor |Ai — A^| there is linear repulsion of complex eigenvalues with 
distance from the real axis, while in the qu-r case there is quadratic repulsion 
from the real axis due to an additional factor |Ai — A^| coming from f q (X\) 2 . 

Correlation functions as Pfaffians It is convenient to abandon the 
Aj's in favour of the real 2-dimensional vectors Zj = Xj with Zj = Xj + iyj 
and the two dimensional volume element d 2 z = dxdy (in contrast, in (|18.5.28|) 
dAdA* = — 2idxdy). Equation (118.5.280 defines the symmetrized jpdf of eigen- 
values P(zi, . . . , zn) in the obvious way (however, care must be taken to account 
for complex conjugate pairs and the ordering of eigenvalues before symmetriz- 
ing the density). P{z\, . . . , z^) is a formal density, as it contains delta-functions 
accounting for complex conjugate pairs and also for eigenvalues on the real line. 

The eigenvalue correlation functions (|18.2.5p (more precisely they are mea- 
sures, since they contain delta-function contributions) can be obtained by func- 
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tional derivatives from the generating functional [Tra98 , KanQH ISom08| 

Z[g] = [ d 2 Zl ... f d 2 z N g{z x ) . . . g{z N ) P(z 1} z 2 , z N ) (18.5.30) 



as R n (z 1 ,z 2 , ...,z n ) = 5g{zi) 5 " Sg(Zn) Z[g] \ g{z) = r They can also be obtained by 
integration from the symmetrized jpdf via (|18.2.5p . This implies a simple in- 
tegration theorem for the full correlations R n (an Pfaffian analogue of Mehta's 
'integrating out' lemma, see (|18.2.9p ). Restricting to the smooth (in the upper 
half-plane) and singular (on the real line) parts of the R n leads to a more sophis- 
ticated integration theorem [Kan05b] connecting these quantities recursively. 
We find it more convenient to work with the generating function Z[g] which can 
be calculated in the form of a Pfaffian with the help of the Grassmann integral 
representation of the Vandermonde determinant |Som08] Tli<i<j<N( z i ~ z j) = 

/ dxi---dxNH i= i(E l k=l z i~ 1 Xk) with XkXi = -XiXk- Since the integrand 
g(z\) . . . g(z^) in (|18.5.30p is symmetric in the Zi's, one can integrate it against 
the joint distribution of ordered eigenvalues as in (|18.5.28|) in order to obtain 
Z[g\. This helps to avoid the absolute values in the Vandermonde determinant, 
making it possible to apply the Grassmann integral above. Recalling the Grass- 
mann integral representation Pfaff (A H ) = J dxi ■ ■ ■ dxN exp(-± Y^ki XkMiXi) 
for the Pfaffian of an antisymmetric matrix (Am) one finds immediately 

Z[g] = Cat Pfaff (A kl ), (18.5.31) 

with A kl = J Jd 2 z 1 d 2 z 2 J : (z 1 , z 2 ) z^zlf 1 g(z 1 )g(z 2 ) and 

F(z\,z 2 ) = f(z 1 )f(z 2 )(2i5 2 (z 1 - z 2 )sgn(yi) + 5(y 1 )5(y 2 ) sgn(x 2 - x x ) ) . 

(18.5.32) 

This simple form of J- is a consequence of the fact that in (|18.5.28p the normal- 
ization constant does not depend on the chosen number m of complex conjugate 
pairs. The first and second terms on the rhs account for eigenvalues coming 
in complex conjugate pairs and eigenvalues on the real axis, respectively. The 
second term vanishes in the qu-r case since f q = on the real axis. 

Putting g(z) = l + u(z) and expanding Z[l + u] in powers of u one again ob- 
tains a series of Pfaffians [Som08j. This method goes back to Mehta's alternate 
variables [Meh04] and Tracy and Widom's paper [Tra98] . With A k i\ g=l = Am, 
A = (Aki) and defining the kernel 

N N 

JC N ( Zl ,z 2 ) =J2J2 A M z t 1 4~ 1 (18.5.33) 

k=l 1=1 

the n-point densities are given by 

R n (z 1 , ...,z n ) = PfafT(Qw) with Q M = ( ^ ^ ) . (18.5.34) 
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This means the Pfaffian of the 2n x 2n matrix built of the n 2 quaternions 
(2x2 matrices) Q k i, k, I = 1, 2, . . . , n with entries G k i = f d 2 z JCiv(z k , z)J r (z, Z[), 
K k i = K N (z k ,zi), and W k i = fd 2 zf d 2 z' T{z k , z)K N {z, z')T{z' , zf) -T(z k ,zi) 
and symmetries K k \ = —Ki k , W k i = —Wi k . With the above expressions one 
finds, e.g., the one-point density: 

Riizt) = J d 2 z 2 T(zi,z 2 )K N (z 2 , z x ) = R?(z) + S(y)R{ H (x) (18.5.35) 

and the two-point density, see Eq. (21) in [Som07]. Note that T{z\,z-i) is 
composed of two parts. Correspondingly, R\(z) contains two parts, a smooth 
part Ri(z) which describes the density of complex eigenvalues and a singular 
part Rf (x) that describes the density of real eigenvalues. In the qu-r case 
Ri(x) = and only the complex part remains. Restricting oneself to points 
z\ , Zi , . . . , z n in the upper half-plane no terms containing delta functions like 
5 2 {zi—z*) appear in the correlation functions and (|18.5,34p gives the correlations 
directly in terms of Pfaffians involving the kernel /Cat- Also, the real correlations 
can be obtained by considering only the terms concentrated on the real axis 
|Som07l IFor07| . These however involve still some real integrations [Som07] . 

Kernel, characteristic polynomials In order to use (|18.5.34j) one needs 
to know the kernel 1Cn(zi, z 2 ) which we now are going to find. It follows from 
(118.5.351) and (I18.5.32|) that R c (z) = 2f(z)f(z*)\IC N (z, z*)\. On the other hand 
one can find R\{z) (and R\{z)) directly from the jpdf ( I18.5.28[) by integrating 
out N — 1 variables. Since the eigenvalues are real or come in complex conjugate 
pairs this leads to Ri(z) oc \z — z*\(det(J — z) det( J — z*))n-2, where (. . -)n-2 
means averaging over the ensemble (|18.5.24[ ) in N — 2 dimensions [Ede97j . On 
comparing the two expressions for Rf (z) one arrives at the important relation!! 

JC N (z,z*) oc (z - z*)(det(J - z) det(J - z*)) N _ 2 ■ (18.5.36) 

Thus, we are left with the task of calculating (det(J — u) det(J-v)) n- Writing 
the determinants as Grassmann integrals, det(J — u) = J TJr]e~ 11 ^ J ~ u ^ 11 , 

P N (u,v) = (det( J - u) det( J - v)) N = J Dr/D( e^***^ e5« Tr J(v^H( r ) 2 )N 

(18.5.37) 

where we have used the Gaussian property. The lowest cumulants of J are 
(JijjN = and (JijJ k i)N = AjfcAjj, where in the real case Ai k = 5%k, while in 
the qu-r case Ai k = Zi k with Z being the symplectic unit. 

2 One can also express the kernel and/or correlation functions via averages of the char- 
acteristic polynomials in the qu-r and complex Ginibre ensembles and in some ensem- 
bles beyond Gaussian. The averages like on the rhs in (|18.5.36[) can be computed in a 
variety of ways and this gives and an alternative way of calculating the eigenvalue cor- 
relation functions in the complex plane in a variety of random matrix ensembles, see 
|Ede94l |Fyo99| IAke03l lAke07l |Fyo07 | IAke09a] , 
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In the real case, after a complex Hubbard-Stratonovich (HS) transformation, 
y(u,v)= e-^ 2 (^J dri*d V dCdCe uv * v+vC ^ ia ^^ +idv ^ . (18.5.38) 



On evaluating the simple Grassmann integration, one arrives at 

™ e -H\ uv + \ a \Y = N\y ^— /- . (18.5.39) 
7r ' n! 

n=0 

After restoring the normalization one finds the kernel 

N ~ 2 I \n 

^E^^ - ,, . ( 18 ,,„, 

n=0 

This gives immediately the density of complex eigenvalues [Ede97] , cf (I18.2.11|) , 
R°(z) = M e 2 ^erfc(^|y|) ^ ~ ± ^ , z = * + h, . (18.5.41) 
and, in view of (|18.5.35 1) . the density of real eigenvalues [Ede94j 

flf (l) = 7SW - i) + 2"-/*r ( iv/2r (~ -t) ■ (18 - 5 ' 42) 

with j*(N,x) = x- N {l-T(N,x)/T(N)). The rhs in (118.5.411) and (|18.5.42l) is 
analytic in TV and these equations hold for odd N as well [Som08, F_or09aJ . 

One can easily analyze (|18.5.4ip and (|18.5.42[) in the limit N — > oo. For 
example, the average number (n^)jv of real eigenvalues can be found by inte- 
gration, recovering the result by Edelmann, Kostlan and Shub |Ede94] 



= i+ ~L n-tW + t) - v - (18 - 5 - 43) 

and verifying the conjecture (nji)N oc y~N made in [Leh91| . Interestingly, the 
variance of tir is also proportional to y/~N [For 07]. For large N and away from 
the real line (|y| 3> 1) the density of complex eigenvalues obeys the circular 
law (|18.2.1ip with the same edge profile (|18.2. 12 j) as in the complex Ginibre 
ensemble. Inside the circle \z\ < yN and close to the real line (y = O(l)) 
Ri(z) ~ y^/V |y| e 2y erfc(\/2|y|). One also finds a constant density of real 
eigenvalues inside the same circle [Ede94] i?^(x) ~ Ty^®^^^ — wmcn is 
consistent with the asymptotics in (|18.5.43p . This constant density is in contrast 
to the Wigner Semicircle Law for Gaussian Hermitian or real symmetric matri- 
ces. In the transitional region around the end point x = y/~N |For07} IBor09| 

R?(VN + it) ~ — — erfc(V2«) H -7= exp(-u 2 ) erfc(-ti) . 

2^ 4^F 
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Although in the qu-r case £#(21,22) can De obtained as saddle point integral 
similarly to (|18.5,39p . it is convenient to derive it directly from the integral Am 

A kl = j d 2 z 2(z - **)e-W V- 1 ** 1 " 1 = 2n(k\ S kil ^ - 11 <5^_i) . (18.5.44) 

We see that in the qu-r case the matrix (Am) has a simple tridiagonal structure, 
while in the real case its inverse is tridiagonal (as evident from (|18.5,40p ). Due 
to this structure we were able to find a beautiful formula for all Schur function 
averages (moments) in the real Ginibre ensemble }Som09j , With this formula 
one can calculate the moments of symmetric functions in eigenvalues by making 
use of the Schur function expansion. A similar formula exists for the qu-r 
Ginibre ensemble |For09b| . 

Using the duplication formula for T(z) and introducing a k = 2 k / 2 T(k/2), 

A k i = — V 71 /^ a fc+i a 2+i e ki in the qu-r case, A k \ = a k ai e H l in the real case, 

where (e k i) = tridiag(l,0, —1) is the antisymmetric tridiagonal matrix with l's 
below the main diagonal and -l's above and its inverse (also antisymmetric) 
e^ 1 = 1 for k odd and j even in the upper triangle k < j and the remaining 
entries being zero in that triangle. From these formulas and (|18,5.33p one finds 
the kernel £#(#1,22) in the qu-r case thus recovering Mehta's result [Meh04] 

I N N 2~( fc+i )/ 2 

Also, the normalization constant Cat in (|18.5.28p follows via (|18.5.31|) 

[f ) N/4 Uk=i afc+i = W' 2 I' 3L.(iV - 1)! in the qu-r case 
fjf =1 a k = (2 v / 27t) jv / 2 0! 2\...{N — 2)! in the real case. 

(18.5.45) 

We see that the problem with all correlations and moments is to calculate 
the matrix A k i and its inverse A7, . In the qu-r case the calculation of A k i 
is simple and leads then to A7, . In the real case the calculation of A7^ via 
characteristic polynomials is simple and a direct calculation of A k \ and taking 
then the inverse is much more involved but possible with the help of the method 
of skew orthogonal polynomials wrt to the form T{z\, Z2), as was demonstrated 
in [For07t IBor09| . The method of skew-orthogonal polynomials can also be used 
to derive the correlation functions in the qu-r ensemble [Kan02| . We will quote 
these results in connection with more general elliptic ensembles. 

18.6 Real and qu-r elliptic ensembles 

Consider now the two families of normalized measures 

dfi T {J) = B T e~^^) Tl{JJf ^ iJ2+Ji2)) |DJ| (18.6.46) 
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on the space oi N x N matrices subject to the symmetry constraints ([18,5,251) . 
The normalization constant B T is (1 — r)^ N<yN ^ l ^ 2 {l -\-T)~ N<yN+l ^ 2 in the real 
case and (1 + t)~ n ^ n ~ 1 ^ 2 {1 — t)~ n ( n+1 ^ 2 in the qu-r case. When r varies 
between —1 and +1 d/u T interpolates between anti-Hermitian (antisymmet- 
ric or antiselfdual, r = —1) and Hermitian (symmetric or selfdual, r = +1) 
ensembles. On comparing to (j!8.5.24|) . one concludes that the jpdf of eigenval- 
ues in ensembles (118. 6. 46ft can simply be obtained by a rescaling of the term 
Tr J jt and otherwise multiplying with the factor exp(A 2 + A* 2 )t/4(1 — r 2 )), 
thus leading to the same expression (118.5.281) with /(A) replaced by / T (A), with 
/ 2 (A) = / 2 (A/\/l — t 2 ) exp{ 2(i- T $) + ^* 2 )} an d a new normalization con- 
stant Cn,t- Correspondingly, the skew-symmetric form is given by (|18.5.32f) . 
again with f(z) replaced by f T (z). Hereby the matrices A^, Aj^ , the kernel 
JCm{z\, z 2 ) and the correlations can be obtained. 

Kernel in the elliptic case The kernel is again related to the average 
of two characteristic polynomials (|18.5.36j) (jl8.5.37f) . but now the second cu- 
mulants of J are {JijJkijN = ^-ik^-jl + T ^U^-jk with Ajfc as before. 

In the real case one can evaluate the corresponding Grassmann integral by 
introducing two complex and two real HS transformations. This gives 

P N (u,v) = [ ^^e-H 2 -H 2 / ^V^ 2 )/ 2 /f T (18.6.47) 
J TT 7T J 2ir 

where P r)T is the Pfaffian of a 4-dimensional antisymmetric matrix equivalent 
to a 4-fold Gaussian Grassmann integral with the value P r>T = (u+ ix- v /r)(w + 
iyv^) + |a| 2 + r|6| 2 . From (118.6.471) one can obtain [Ake09 a] the kernel 

JC N ( Zl ,z 2 ) = V p^z^^L-P^zz) (18.6.48) 

k,l=l 

expanded in terms of skew orthogonal monic polynomials Pk~\{z) defined by 

J d 2 Zl d 2 z 2 F(z 1 , z 2 )Pk-i(zi)Pi- 1 (z 2 ) = Z u r,_i = r fe _! Z kl . (18.6.49) 

For r > 0, (|18.6.48[) was also obtained directly from the jpdf in |For08| . In 
this case the skew-orthogonal polynomials can be expressed [For08] in terms of 
the Hermite polynomials P 2n (z) = p 2n (z) and P 2n +i = P2n+i(z) - 2np 2n ^i(z) 
with Pn(z) = (t /2) n l 2 H n (z I 'y/2r) being the scaled Hermite polynomials as in 
Section ["18.41 and with r n = r n+ i = 2^/2irn\ (1 + r) for n even. 

In the qu-r case the saddle point integral representation of the kernel fol- 
lowing from its relation to the average of characteristic polynomials is more 
involved, however the kernel can again be expanded in skew orthogonal polyno- 
mials (for t > 0) [Kan02]. The resulting expression is the same as in (|18.6.48[) 
only now P 2n (z) = Yd=o ^qrP2l(z) and P 2n+ i = p 2n+ i(z), with the polynomials 
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p n as before and the constants r n in (|18.6.49p given by r n = r n _i = 2tt n\ (1— r) 
for n odd. Note that the expressions for r n give immediately the normaliza- 
tion constant in (|18.5.28|) : Cm,t = (1 + t)~ n / 2 Cn in the real case and and 
Cn,t = (1 - t)~ n/2 C n in the qu-r case with C N as in (|18.5.45p . 

Having an explicit form for the kernel, all correlations, including in the real 
case the real-real and real-complex ones, can be found I 'orOsl lBor091 ISom 08j . 

Strongly non-Hermitian limit In this limit u, v are assumed to be of 
order ^/N as N — > oo. We are going to evaluate the kernel )Cn(u,v) (|18.5.36|) 
by the saddle point analysis of Pn-2(u,v). Let us start with the circular real 
case (|18.5.39p which is the simplest. Then the saddle point equations are a* = 
Na , and a = — ^ a , with two solutions a = a* = and aa* = N — uv. 

uv+aa* uv+aa* 

The first saddle-point yields Pn-2 — (uv) N ~ 2 /{\ — N/uv) and the second, 
which is actually a manifold, yields Pn-2 — N N ~ 2 \/2irN e uv ~ N . Thus the 
second is dominating for Re(m>) < N and otherwise the first is (excluding the 
neighborhood of N = uv). Hence, the 1-point density R\{z) is asymptotically 
R\(z) = l/ir for \z\ < \^N and is exponentially small outside this circle. Also, 
this gives immediately the asymptotic form of the kernel in the circular case in 
the bulk (Re0i2 2 ) < N): 

JC N (zi,z 2 ) ~ (2v / 2^) _1 (zx - z 2 )e zlZ2 as N oo. (18.6.50) 

Obviously, (|18.6.50p could have been obtained directly from (I18.5.40p . Here 
we have determined in addition the region of validity. The correlations (re- 
point densities) are then given with this formula in the region \zj\ < yN , 
j = 1, 2, . . . , re. Hence we call it the circular real case. 
By similar reasoning we find in the circular qu-r case 

KL N {z l ,z 2 )~{2Tr{z2-z 1 ))- 1 e ZlZ2 as iV -> oo, (18.6.51) 

for Re(ziZ2) < N with both Z\ and z 2 (and also Z\ — z 2 ) being of order y/~N. 
This implies the circular law in the qu-r case: R\{z) ~ 1/tt inside the circle 
\z\ = y/~N and away from the real line (Imz ~ y/~N) and R\{z) ~ outside. 

For the elliptic real case the relevant saddle point leads to the condition for 
the elliptic support of R\{z). This has been found numerically and analytically 
using the replica trick already in the early paper |Som88] . Similarly, for the 
elliptic qu-r case we expect it can be shown by a saddle point analysis that 
in the large N- limit R\{z) is constant inside the ellipse with main half-axes 
VN(1 ± t) along the x-and y- direction, and zero outside. 

Having the expressions (|18.6.48[ ) for the kernel in terms of skew orthogonal 
polynomials, one can find the scaling limit of the kernel and the correlations for 
iV to infinity considering z\ 2 as being of order 1 and letting the edge going to 
infinity. Since the 1-point density is constant it is not necessary to unfold. It 
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turns out that in the real case for r = this limit CqoO^i, £2) is already given 
by expression (|18.6.50p . For r > K.qo(zi, &z) has been calculated in [For08j 



Z\ — Z% ( Z\Z2 T[z\ + Z%) 



with the normalization corrected. Obviously this is also valid for —1 < r < 0. 
In the circular qu-r case (r = 0) the asymptotic form of the kernel has been 
calculated in [Kan02] and is already found in [MehQ4j in a somewhat different 
form. Scaling it appropriately one obtains for the elliptic qu-r case 

K^j^—J exp ( S-^ \ erf ( ~ " 2 =) . (18.6.53) 

Both expressions (|18.6.52l [15X33"]) yield the bulk density Ri{z) ~ l/vr(l - r 2 ) 
corresponding to an elliptic support with main half axes y/N (l±r). In the real 
case the density of real eigenvalues is asymptotically again constant Rf(x) ~ 
1/a/2"7t(1 — r 2 ) which together with the support |x| < y/~N (1+r) gives a number 
of real eigenvalues of ^2^(1 + r)/vr(l - r) (cf (|18.5.43|) 1 |For08j . 

Weakly non-Hermitian limit Finally we consider the limits of weak 
non-Hermiticity |Fyo97| . Here we put r = 1 — a 2 with 1 — r being of order 1/N 
and let N go to infinity assuming Zi to be in the neighborhood of the origin. In 
the real case in this limit the kernel is given by [For08j 



N f 1 

Kn(zi,Z2) — — / du u ex.p(—Na 2 u 2 )sm(VNu(zi — Z2)), (18.6.54) 
2vr Jo 

and f 2 (z) ~ erfc( ^2a ' )• I n the- qu-r case f 2 {z) ~ e 5^ and |Kan02| 

K N (zi,Z2) ~ - / — exp(-ATaV) si^x/iV ufo - z 2 )) . (18.6.55) 

Aa 6 N 6 / z Jo u 

As in the complex elliptic ensemble, it is convenient to scale Z\,z%,a and /Cat 
with the local mean level spacing (here = tt/^/N since the ensembles go in 
this limit to GOE/GSE with the semicircular density of eigenvalues p sc (x) = 
l^/N -x 2 /A). Then one obtains a universal form of the correlations. For 
example, on the large rc-scale (z = x + ly) the 1-point density can be written as 
R\{z) ~ p sc (x) 2 P(p sc (x)y, p sc (x)a) with the probability density function 

P(y,a) = 8{y) j\u e -^ W + vrerfc £ du u e -^ W sinh(2^|y|) 

in the real case (GOE limit) [Efe97] and in the qu-r case (GSE limit) [Kol99j 

P(y,a)= J exp ( -^j ) f — exp(— w 2 a 2 u 2 ) s'mh(2iruy) . 
ir- > ' z a' > \ a J Jo u 
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Note that in the qu-r case despite the fact that the complex eigenvalue den- 
sity vanishes exactly on the real axis in the limit pairs of complex conjugate 
eigenvalues approach the real axis and collapse giving Kramers degeneracy for 
the Hermitian ensemble with symplectic symmetry. In the real case despite the 
fact that the density of real eigenvalues becomes constant, but very low, com- 
plex eigenvalues approach the real axis and in the limit give rise to the Wigner 
semicircle density which can be considered as a projection of the elliptic law in 
the complex plane onto the real axis. 
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